Analytical validation and performance characteristics of a 48-gene next-generation sequencing panel for detecting potentially actionable genomic alterations in myeloid neoplasms.

Identification of genomic mutations by molecular testing plays an important role in diagnosis, prognosis, and treatment of myeloid neoplasms. Next-generation sequencing (NGS) is an efficient method for simultaneous detection of clinically significant genomic mutations with high sensitivity. Various NGS based in-house developed and commercial myeloid neoplasm panels have been integrated into routine clinical practice. However, some genes frequently mutated in myeloid malignancies are particularly difficult to sequence with NGS panels (e.g., CEBPA, CARL, and FLT3). We report development and validation of a 48-gene NGS panel that includes genes that are technically challenging for molecular profiling of myeloid neoplasms including acute myeloid leukemia (AML), myelodysplastic syndrome (MDS), and myeloproliferative neoplasms (MPN). Target regions were captured by hybridization with complementary biotinylated DNA baits, and NGS was performed on an Illumina NextSeq500 instrument. A bioinformatics pipeline that was developed in-house was used to detect single nucleotide variations (SNVs), insertions/deletions (indels), and FLT3 internal tandem duplications (FLT3-ITD). An analytical validation study was performed on 184 unique specimens for variants with allele frequencies ≥5%. Variants identified by the 48-gene panel were compared to those identified by a 35-gene hematologic neoplasms panel using an additional 137 unique specimens. The developed assay was applied to a large cohort (n = 2,053) of patients with suspected myeloid neoplasms. Analytical validation yielded 99.6% sensitivity (95% CI: 98.9-99.9%) and 100% specificity (95% CI: 100%). Concordance of variants detected by the 2 tested panels was 100%. Among patients with suspected myeloid neoplasms (n = 2,053), 54.5% patients harbored at least one clinically significant mutation: 77% in AML patients, 48% in MDS, and 45% in MPN. Together, these findings demonstrate that the assay can identify mutations associated with diagnosis, prognosis, and treatment options of myeloid neoplasms even in technically challenging genes.

Genome sequencing of 18 francisella strains to aid in assay development and testing.

Francisella tularensis is a highly infectious bacterium with the potential to cause high fatality rates if infections are untreated. To aid in the development of rapid and accurate detection assays, we have sequenced and annotated the genomes of 18 F. tularensis and Francisella philomiragia strains.

Accurate read-based metagenome characterization using a hierarchical suite of unique signatures.

A major challenge in the field of shotgun metagenomics is the accurate identification of organisms present within a microbial community, based on classification of short sequence reads. Though existing microbial community profiling methods have attempted to rapidly classify the millions of reads output from modern sequencers, the combination of incomplete databases, similarity among otherwise divergent genomes, errors and biases in sequencing technologies, and the large volumes of sequencing data required for metagenome sequencing has led to unacceptably high false discovery rates (FDR). Here, we present the application of a novel, gene-independent and signature-based metagenomic taxonomic profiling method with significantly and consistently smaller FDR than any other available method. Our algorithm circumvents false positives using a series of non-redundant signature databases and examines Genomic Origins Through Taxonomic CHAllenge (GOTTCHA). GOTTCHA was tested and validated on 20 synthetic and mock datasets ranging in community composition and complexity, was applied successfully to data generated from spiked environmental and clinical samples, and robustly demonstrates superior performance compared with other available tools.

Draft genome sequence of the rubber tree Hevea brasiliensis.

BACKGROUND: Hevea brasiliensis, a member of the Euphorbiaceae family, is the major commercial source of natural rubber (NR). NR is a latex polymer with high elasticity, flexibility, and resilience that has played a critical role in the world economy since 1876. RESULTS: Here, we report the draft genome sequence of H. brasiliensis. The assembly spans ~1.1 Gb of the estimated 2.15 Gb haploid genome. Overall, ~78% of the genome was identified as repetitive DNA. Gene prediction shows 68,955 gene models, of which 12.7% are unique to Hevea. Most of the key genes associated with rubber biosynthesis, rubberwood formation, disease resistance, and allergenicity have been identified. CONCLUSIONS: The knowledge gained from this genome sequence will aid in the future development of high-yielding clones to keep up with the ever increasing need for natural rubber.

Bacterial genome annotation.

Annotation of prokaryotic sequences can be separated into structural and functional annotation. Structural annotation is dependent on algorithmic interrogation of experimental evidence to discover the physical characteristics of a gene. This is done in an effort to construct accurate gene models, so understanding function or evolution of genes among organisms is not impeded. Functional annotation is dependent on sequence similarity to other known genes or proteins in an effort to assess the function of the gene. Combining structural and functional annotation across genomes in a comparative manner promotes higher levels of accurate annotation as well as an advanced understanding of genome evolution. As the availability of bacterial sequences increases and annotation methods improve, the value of comparative annotation will increase.

Genome of the cyanobacterium Microcoleus vaginatus FGP-2, a photosynthetic ecosystem engineer of arid land soil biocrusts worldwide.

The filamentous cyanobacterium Microcoleus vaginatusis found in arid land soils worldwide. The genome of M. vaginatus strain FGP-2 allows exploration of genes involved in photosynthesis, desiccation tolerance, alkane production, and other features contributing to this organism's ability to function as a major component of biological soil crusts in arid lands.

Genomics for key players in the N cycle from guinea pigs to the next frontier.

While sequencing methods were available in the late 1970s, it was not until the human genome project and a significant influx of funds for such research that this technology became high throughput. The fields of microbiology and microbial ecology, among many others, have been tremendously impacted over the years, to such an extent that the determination of complete microbial genome sequences is now commonplace. Given the lower costs of next-generation sequencing platforms, even small laboratories from around the world will be able to generate millions of base pairs of data, equivalent to entire genomes worth of sequence information. With this prospect just around the corner, it is timely to provide an overview of the genomics process: from sample preparation to some of the analytical methods used to gain functional knowledge from sequence information.

Tessellation analysis of glomerular spatial arrangement in mice with heritable renal hypoplasia.

Renal hypoplasia results from an insufficient kidney volume caused, in part, by a deficient number of glomeruli. The purpose of this study was to apply tessellation analysis to determine whether glomerular point patterns differed between adult normal (WT) and mutant (Br) mice with heritable renal hypoplasia and to delineate a spatial distribution accounting for the observed patterns. Kidneys from adult WT and Br mice were collected, processed with routine light histology and representative transverse sections were photographed. Cortical area and perimeter were calculated from traced tissue contours and glomeruli were identified and digitized. Voronoi tessellations were constructed and average parameters for Voronoi polygon number, area, perimeter and edge counts as well as spatial metrics comprising nearest neighbor and centroidal distances were calculated and compared. Point distributions were simulated by randomizing glomerular coordinates from each section and plotting the new points utilizing uniform random, Gaussian random, or isotropic functions. Average nearest neighbor distances were generated for each specimen and ranked with respect to corresponding values generated from 1,000 iterations for each simulated set. Results showed that WT and Br were significantly different for each parameter suggesting that WT kidneys possessed more glomeruli, but these were less clustered compared to Br. Simulations suggested that WT and Br demonstrated similar, but not identical, underlying glomerular spatial distributions. Defective gene expression in Br is important for determining glomerular number and the defective pattern likely results from a heterochronic disturbance consisting of a truncated growth trajectory during embryonic kidney development.

Globins synthesize the second messenger bis-(3'-5')-cyclic diguanosine monophosphate in bacteria.

Globin-coupled sensors are heme-binding signal transducers in Bacteria and Archaea in which an N-terminal globin controls the activity of a variable C-terminal domain. Here, we report that BpeGReg, a globin-coupled diguanylate cyclase from the whooping cough pathogen Bordetella pertussis, synthesizes the second messenger bis-(3'-5')-cyclic diguanosine monophosphate (c-di-GMP) upon oxygen binding. Expression of BpeGReg in Salmonella typhimurium enhances biofilm formation, while knockout of the BpeGReg gene of B. pertussis results in decreased biofilm formation. These results represent the first identification a signal ligand for any diguanylate cyclase and provide definitive experimental evidence that a globin-coupled sensor regulates c-di-GMP synthesis and biofilm formation. We propose that the synthesis of c-di-GMP by globin sensors is a widespread phenomenon in bacteria.

Cloning, expression, and purification of the N-terminal heme-binding domain of globin-coupled sensors.

Globin-coupled sensors (GCSs) are multidomain proteins, consisting of an N-terminal globin domain fused to a variety of C-terminal transmitter domains. Functional classification of GCSs is based on the transmitter domain(s) they possess, broadly falling under either aerotaxis or gene regulation. This chapter describes methods and strategies for cloning, expression, and purification of GCSs for spectroscopic analysis and determination of the minimum heme-binding domain, together with bioinformatic approaches for database searching and examination of domain architectures.

Molecular dynamics and in vitro analysis of Connexin43: A new 14-3-3 mode-1 interacting protein.

The interaction of cellular proteins with the gap junction protein Connexin43 (Cx43) is thought to form a dynamic scaffolding complex that functions as a platform for the assembly of signaling, structural, and cytoskeletal proteins. A high stringency Scansite search of rat Cx43 identified the motif containing Ser373 (S373) as a 14-3-3 binding site. The S373 motif and the second best mode-1 motif, containing Ser244 (S244), are conserved in rat, mouse, human, chicken, and bovine, but not in Xenopus or zebrafish Cx43. Docking studies of a mouse/rat 14-3-3 homology model with the modeled phosphorylated S373 or S244 peptide ligands or their serine-to-alanine mutants, S373A or S244A, revealed that the pS373 motif facilitated a greater number of intermolecular contacts than the pS244 motif, thus supporting a stronger 14-3-3 binding interaction with the pS373 motif. The alanine substitution also reduced more than half the number of intermolecular contacts between 14-3-3 and the S373 motif, emphasizing the phosphorylation dependence of this interaction. Furthermore, the ability of the wild-type or the S244A GST-Cx43 C-terminal fusion protein, but not the S373A fusion protein, to interact with either 14-3-3 or 14-3-3zeta in GST pull-down experiments clearly demonstrated that the S373 motif mediates the direct interaction between Cx43 and 14-3-3 proteins. Blocking growth factor-induced Akt activation and presumably any Akt-mediated phosphorylation of the S373 motif in ROSE 199 cells did not prevent the down-regulation of Cx43-mediated cell-cell communication, suggesting that an Akt-mediated interaction with 14-3-3 was not involved in the disruption of Cx43 function.

Globin-coupled sensors, protoglobins, and the last universal common ancestor.

The strategy for detecting oxygen, carbon monoxide, nitric oxide, and sulfides is predominantly through heme-based sensors utilizing either a globin domain or a PAS domain. Whereas PAS domains bind various cofactors, globins bind only heme. Globin-coupled sensors (GCSs) were first described as regulators of the aerotactic responses in Bacillus subtilis and Halobacterium salinarum. GCSs were also identified in diverse microorganisms that appear to have roles in regulating gene expression. Functional and evolutionary analyses of the GCSs, their protoglobin ancestor, and their relationship to the last universal common ancestor (LUCA) are discussed in the context of globin-based signal transduction.

Genome sequence of the deep-sea gamma-proteobacterium Idiomarina loihiensis reveals amino acid fermentation as a source of carbon and energy.

We report the complete genome sequence of the deep-sea gamma-proteobacterium, Idiomarina loihiensis, isolated recently from a hydrothermal vent at 1,300-m depth on the Loihi submarine volcano, Hawaii. The I. loihiensis genome comprises a single chromosome of 2,839,318 base pairs, encoding 2,640 proteins, four rRNA operons, and 56 tRNA genes. A comparison of I. loihiensis to the genomes of other gamma-proteobacteria reveals abundance of amino acid transport and degradation enzymes, but a loss of sugar transport systems and certain enzymes of sugar metabolism. This finding suggests that I. loihiensis relies primarily on amino acid catabolism, rather than on sugar fermentation, for carbon and energy. Enzymes for biosynthesis of purines, pyrimidines, the majority of amino acids, and coenzymes are encoded in the genome, but biosynthetic pathways for Leu, Ile, Val, Thr, and Met are incomplete. Auxotrophy for Val and Thr was confirmed by in vivo experiments. The I. loihiensis genome contains a cluster of 32 genes encoding enzymes for exopolysaccharide and capsular polysaccharide synthesis. It also encodes diverse peptidases, a variety of peptide and amino acid uptake systems, and versatile signal transduction machinery. We propose that the source of amino acids for I. loihiensis growth are the proteinaceous particles present in the deep sea hydrothermal vent waters. I. loihiensis would colonize these particles by using the secreted exopolysaccharide, digest these proteins, and metabolize the resulting peptides and amino acids. In summary, the I. loihiensis genome reveals an integrated mechanism of metabolic adaptation to the constantly changing deep-sea hydrothermal ecosystem.

Ancestral hemoglobins in Archaea.

Hemoglobins are ubiquitous in Eukarya and Bacteria but, until now, have not been found in Archaea. A phylogenetic analysis of the recently revealed microbial family of globin-coupled heme-based sensors suggests that these sensors descended from an ancient globin-only progenitor, or a protoglobin (Pgb). Here, we report the discovery and characterization of two Pgbs from the Archaea: ApPgb from the obligately aerobic hyperthermophile Aeropyrum pernix, and MaPgb from the strictly anaerobic methanogen Methanosarcina acetivorans. Both ApPgb and MaPgb bind molecular oxygen, nitric oxide, and carbon monoxide by means of a heme moiety that is coordinated to the protein through the F8 histidine (histidine 120). We postulate that these archaeal globins are the ancestors of contemporary hemoglobins.

The diversity of globin-coupled sensors.

The recently discovered globin-coupled sensors (GCSs) are heme-containing two-domain transducers distinct from the PAS domain superfamily. We have identified an additional 22 GCSs with varying multi-domain C-terminal transmitters through a search of the complete and incomplete microbial genome datasets. The GCS superfamily is composed of two major subfamilies: the aerotactic and gene regulators. We postulate the existence of protoglobin in Archaea as the predecessor to the chimeric GCS.